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Abstract 


Exploitation  of  the  backscattered  field  polarization  over  the  wide  electromagnetic  spectrum,  from 
visible  to  microwave  frequencies,  provides  an  approach  to  advanced  target  recognition.  It  is 
analytically  proven  that  the  current,  established  framework  for  full  coherent  characterization  of 
the  scattered  field  maximizes  the  extracted  target  information.  It  is  also  shown  that  such  a 
methodology,  which  is  theoretically  similar  to  the  concept  of  “partial  or  compact  polarimetry”, 
yields  comparable  results  to  full  or  quadrature-polarized  systems  by  incoiporating  judicious 
assumptions  and  assuming/implementing  optimal  transmitted  or  illumination  field  polarizations. 
On  this  basis,  common  characteristic  features,  interworking  and  fusion  of  different  polarimetric 
sensor  products  in  different  regions  of  spectrum,  e.g.,  radar/S AR  and  Electro-Optical,  are 
investigated  and  formulated  within  a  robust  framework  based  on  full  coherent  characterization  of 
the  scattered  field. 


Resume 


L’exploitation  de  la  polarisation  de  champ  retrodiffusee  sur  une  vaste  plage  du  spectre 
electromagnetique  allant  de  la  region  visible  jusqu’aux  frequences  micro-ondes  donne  une 
approche  de  la  reconnaissance  evoluee  des  objectifs.  On  a  demontre  analytiquement  que  le  cadre 
etabli  actuel  de  caracterisation  coherente  integrale  du  champ  diffuse  permet  d’extraire  un 
maximum  d’ informations  sur  l’objectif.  Cette  demarche  montre  aussi  qu’une  telle  methode  qui  est 
theoriquement  similaire  au  concept  de  la  polarimetrie  partielle  dite  compacte  donne  des  resultats 
comparables  aux  systemes  a  polarimetrie  complete  ou  en  quadrature  en  integrant  de  judicieuses 
hypotheses  et  en  presumant  ou  mettant  en  oeuvre  la  polarisation  optimale  des  champs  transmis  ou 
leur  illumination.  En  se  basant  sur  ce  qui  precede,  on  etudie  les  caracteristiques  communes, 
l’interfonctionnement  et  la  fusion  de  divers  produits  de  capteurs  polarimetriques  dans  diverses 
regions  du  spectre,  p.  ex.  le  radar  classique  ou  le  radar  a  synthese  d’ouverture  et  des  capteurs 
electro-optiques,  et  on  formule  des  hypotheses  pertinentes  dans  un  cadre  robuste  base  sur  la 
caracterisation  coherente  integrale  du  champ  diffuse. 
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Executive  summary 


Polarimetric  Intensity  Parameterization  of  radar  and  other 
remote  sensing  sources  for  advanced  exploitation  and  data 
fusion:  Theory 

R.  Sabry,  P.W.  Vachon;  DRDC  Ottawa  TM  2008-217;  Defence  R&D  Canada  - 
Ottawa;  October  2008. 

Introduction:  Backscattered  Electromagnetic  (EM)  waveform  polarimetry  provides  additional 
information  for  target  recognition  compared  to  conventional  single  channel  remote  sensing 
sources.  The  objective  here  is  to  examine  the  interworking  and  fusion  of  EM  polarimetric  return 
quantities  and  measurables  associated  with  certain  classes  of  targets  in  different  regions  of  EM 
spectrum.  One  prominent  example  is  the  combination  of  radar/SAR  polarimetry  with  EM 
polarimetric  sensors  in  the  electro-optical  domain. 

Results:  Preliminary  study,  formulations  and  simulations  indicate  the  potential  for  advanced 
applications.  Simulated  results  are  obtained  in  the  radar/SAR  domain.  A  general  analytic 
framework,  valid  for  different  regions  of  the  EM  spectrum,  is  established.  Core  conceptual 
similarity  to  “Compact  Polarimetry”  that  is  new,  fast  growing  and  proven  potent  for  advanced 
radar  data  exploitations,  also  indicates  the  significance  of  the  proposed  approach.  Considering  the 
general  property  of  the  scattered  EM  wave  being  polychromatic  and  partially  polarized,  the 
analytic  polarimetric  EM  signal  is  characterized  by  virtue  of  the  Stokes  vector  and  related 
parameters. 

Significance:  The  general  backscattered  polarimetric  wave  in  different  regions  of  the  EM 
spectrum  is  polychromatic  and  partially  polarized.  Also,  it  has  been  proven  that  the  measurement 
potential  of  polarimetric  SAR  in  response  to  a  random  backscatterer  is  maximized  only  if  the  data 
products  are  the  4-vector  Stokes  parameters.  Thus,  using  the  developed  framework  based  on  the 
comprehensive  Stokes  parameters  and  products  provides  a  means  for  robust  sensor  data  fusion. 
The  latter  offers  significant  advancement  in  data  exploitation  and  target  recognition  due  to  the 
following:  target  scattering  characteristics  (i.e.,  identity)  are  better  established  over  a  wider  range 
of  the  EM  spectrum;  and  polarimetric  return  characteristics  for  certain  target  classes  of  interest 
are  consistent  throughout  various  regions  of  the  EM  spectrum  (e.g.,  polarization  selectivity  for 
manmade  objects).  Elence,  proper  fusion  of  EM  polarimetric  returns  in  different  bands  can,  in 
principle,  enhance  the  contrast  of  such  objects/targets  compared  to  those  with  low  or  without  the 
addressed  polarization  selectivity  (e.g.,  clutter).  The  described  EM  polarimetric  nature 
(polarization  selectivity)  manifested  in  different  regions  of  the  spectrum  can  be  effectively 
exploited  for  classification  purposes. 
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Polarimetric  Intensity  Parameterization  of  radar  and  other 
remote  sensing  sources  for  advanced  exploitation  and  data 
fusion:  Theory 

R.  Sabry,  P.W.  Vachon;  DRDC  Ottawa  TM  2008-217;  R  &  D  pour  la  defense 
Canada  -  Ottawa;  Octobre  2008. 

Introduction  :  La  polarimetrie  des  formes  d’ondes  electromagnetiques  retrodiffusees  foumit 
davantage  de  donnees  pour  la  reconnaissance  des  objectifs  que  les  capteurs  de  teledetection 
monovoie  classiques.  Nous  voulons  examiner  ici  l’interfonctionnement  et  la  fusion  des  grandeurs 
et  variables  associees  aux  echos  electromagnetiques  et  polarimetriques  de  certaines  categories 
d’objectifs  dans  diverses  regions  du  spectre  electromagnetique.  Un  exemple  eminent  en  cette 
matiere  est  la  combinaison  de  la  polarimetrie  du  radar  classique  et/ou  du  radar  a  synthese 
d’ouverture  et  des  capteurs  electromagnetiques  polarimetriques  dans  le  domaine  de  l’electro- 
optique. 

Resultats:  Des  etudes,  formulations  et  simulations  preliminaries  montrent  qu’il  existe  un 
potentiel  d’ applications  evoluees.  On  obtient  des  resultats  simules  dans  le  domaine  du  radar 
classique  et/ou  du  radar  a  synthese  d’ouverture.  Un  cadre  analytique  general  valide  pour  diverses 
regions  du  spectre  electromagnetique  est  defini.  11  existe  une  similarite  conceptuelle  fondamentale 
avec  la  «  polarimetrie  compacte  »,  autrement  dit  un  potentiel  d’innovation  eprouve  et  croissant 
rapidement  pour  des  exploitations  evoluees  des  donnees  radar,  ce  qui  demontre  aussi  la  portee  de 
la  demarche  envisagee.  Eu  egard  aux  proprietes  generates  de  l’onde  electromagnetique 
retrodiffusee,  laquelle  est  a  la  fois  polychromatique  et  partiellement  polarisee,  le  signal 
polarimetrique  electromagnetique  analyse  est  caracterise  en  vertu  du  theoreme  de  Stokes  et  des 
parametres  connexes. 

Portee:  L’onde  polarimetrique  retrodiffusee  generale  est  polychromatique  et  partiellement 
polarisee.  On  a  aussi  demontre  que  le  potentiel  de  mesure  du  radar  a  synthese  d’ouverture 
repondant  a  une  retrodiffusion  aleatoire  est  porte  a  son  maximum  uniquement  si  les  donnees 
correspondent  aux  parametres  quadrivectoriels  du  theoreme  de  Stokes.  Done,  l’utilisation  du 
cadre  mis  au  point  fonde  sur  les  parametres  et  les  produits  integraux  du  theoreme  de  Stokes  donne 
un  moyen  de  realiser  une  solide  fusion  des  donnees  de  capteurs.  Cette  fusion  presente  un 
avancement  important  en  matiere  d’ exploitation  des  donnees  et  de  reconnaissance  des  objectifs, 
grace  aux  elements  suivants  :  les  caracteristiques  de  diffusion  de  l’objectif  (c.-a-d.  son  identite) 
sont  mieux  defmies  sur  une  plage  plus  vaste  du  spectre  electromagnetique;  les  caracteristiques 
polarimetriques  des  echos  de  certaines  categories  d’objectifs  qui  nous  interessent  sont  uniformes 
dans  diverses  regions  du  spectre  electromagnetique  (p.  ex.  la  selectivite  de  polarisation  pour  les 
objets  artificiels).  Consequemment,  une  fusion  adequate  des  echos  electromagnetiques 
polarimetriques  dans  diverses  bandes  peut,  en  principe,  ameliorer  le  contraste  de  tels  objets  et/ou 
objectifs,  par  rapport  a  ceux  ou  Ton  a  peu  ou  pas  tenu  compte  de  la  selectivite  de  polarisation  (c.- 
a-d.  du  fouillis  radar).  La  nature  electromagnetique  et  polarimetrique  (de  l’echo)  decrite  plus  haut 
(selectivite  de  polarisation)  qui  se  manifeste  dans  diverses  regions  du  spectre  peut  effectivement 
etre  exploitee  a  des  fins  de  classification. 
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1  Introduction 


Polarization  characteristics  of  the  electromagnetic  (EM)  waveform  scattered  as  a  result  of 
interaction  with  different  objects  and/or  targets  provides  significant  insight  into  the  nature  of 
those  targets,  hence,  providing  valuable  additional  information  for  target  recognition  compared  to 
conventional  remote  sensing  sources.  This  is  supported  by  the  well-established  physics  of  EM 
fields  and  waves  [1  ]-[5] .  Accordingly,  EM  polarimetry,  from  visible  to  microwave  frequencies, 
provides  an  enhanced  capability  for  advanced  target  recognition  in  remote  sensing  applications. 

The  goal  here  is  to  examine  the  interworking  and  fusion  of  EM  polarimetric  observables  and 
measurables  associated  with  certain  classes  of  targets  in  different  regions  of  the  EM  spectrum. 
One  prominent  example  is  the  combination  of  radar/SAR  polarimetry  with  EM  polarimetric 
sensors  in  higher  frequency  bands,  e.g.,  the  visible  or  optical  region.  Characteristic  advantages  of 
such  fusion  or  hybrid  approach  can  be  expressed  as  follows: 

•  First  and  foremost,  in  theory,  target  identity  established  through  scattering 
characteristics  will  be  more  uniquely  defined  in  a  wider  range  of  the  EM  spectrum 
even  if  the  spectrum  expansion  is  incoherent  (i.e.,  incoherent  bandwidths).  Such 
spectral  diversity  provides  a  more  comprehensive  nature  of  scattering  and  the  target. 

•  Polarimetric  backscatter  characteristics  for  a  certain  target  class  (or  classes)  of 

interest  are  consistent  throughout  various  regions  of  the  EM  spectrum  (e.g., 

polarization  selectivity  for  manmade  objects).  Therefore,  proper  fusion  of  EM 
polarimetric  returns  in  different  bands  can,  in  principle,  enhance  the  contrast  of  such 
objects/targets  compared  to  those  with  low  or  without  the  addressed  polarization 
selectivity  (e.g.,  clutter). 

•  The  described  EM  polarimetric  nature  (polarization  selectivity)  manifested  in 

different  regions  of  the  spectrum  can  be  effectively  exploited  for  calibration 
purposes.  For  instance,  specification  of  the  target  polarization  plane  (e.g.,  the 
polarization  orientation  angle)  may  be  extracted  thus  exploiting  data  from  the  sensor 
in  one  domain  for  exploitation  and  calibration  in  the  other  domain.  Such  interworking 
is  applicable  and  effective  for  any  polarization-related  property  associated  with  target 
physics  that  is  consistent  throughout  the  EM  spectrum  (or  essentially  within  the 
operating  bands  of  interest). 

The  described  fusion  of  polarimetric  data  from  SAR  and  other  sensors  from  a  different  region  of 
the  EM  spectrum,  e.g.,  electro-optics  (EO),  should  be  performed  within  a  robust  framework  based 
on  careful  modeling  of  the  EM  wave  scattering  and  characteristic  observables.  As  will  be 
described,  considering  the  general  property  of  the  scattered  EM  wave  being  polychromatic  and 
partially  polarized,  the  analytic  polarimetric  EM  signal  can  be  best-characterized  by  virtue  of  the 
Stokes  vector  and  parameters.  The  scattering  phenomenon  can  also  be  formulated  within  the  same 
framework  using  Mueller,  coherency  or  variance  matrices  [l],[6]-[8]. 
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2  Electromagnetic  (EM)  scattering 


The  identity  of  objects/targets  interacting  with  the  Electromagnetic  (EM)  field  in  remote  sensing 
is  characterized  by  virtue  of  the  EM  scattering  concept.  The  EM  field  scattering  operator  is 
analogous  to  the  conventional  target  reflectivity  function  but  provides  a  broader  and  more 
comprehensive  picture  of  scattering.  In  particular,  it  contains  information  on  the  scattered  field 
polarization  that  is  valuable  for  target  characterization.  Accordingly,  the  polarimetric  scattering 
operator  is  a  Dyadic  or  Tensor  by  nature. 

In  any  region  of  the  EM  spectrum,  target  scattering  can  be  formulated  as: 


Fs(x,  u>)  =  S(x,  u>) .  F;(x,  u>)  (1) 


where  S(x,u;)  denotes  the  scattering  operator  in  the  space  and  time-frequency  domain,  and 
F;.(x,  uj),  Fs(x,  u)  represent  incident  and  scattered  fields,  respectively.  For  typical  applications, 
these  vector  field  components  represent  the  electric  field,  i.e.,  E(  (x,  u)  andE  y(x,u) .  ln(l), 
denotes  the  dot  or  inner  product. 

As  is  evident  from  ( 1 ),  the  identity  of  a  scatterer  at  position  x  is  completely  and  uniquely  defined 
(theoretically)  over  the  whole  spectrum.  Thus,  waveform  diversity  (for  both  incident  and 
scattered  fields)  is  required  to  establish  the  complete  scatterer  identity.  That  is,  target  recognition 
will  be  enhanced  through  spectrum  diversity. 

One  should  note  that  in  the  present  formalism,  the  excitation  (i.e.,  incident  field)  and  receive  (i.e., 
scattered  field)  points  may  be  in  different  space  locations  such  as  those  for  bistatic  radars.  Also, 
the  incident  field  source  of  scattering  may  be  not  controlled  or  coherent. 

The  general  formalism  of  (1),  although  simple  in  conceptual  representation,  is  very  involved  for 
practical  target  recognition  applications.  The  coherent  scattering  target  operator  or  matrix  that 
contains  comprehensive  scatterer  information  may  be  difficult  to  relate  to  target  characteristics 
for  practical  applications.  This  is  due  to  many  factors  such  as  the  noted  incomplete  spectral 
content,  random  effects,  partial  or  non  coherence  (with  respect  to  the  source  excitation/incident 
field)  associated  with  the  scattered  field  or  wave.  Considering  the  above,  the  scattered  field  or 
wave  analysis  becomes  a  viable  choice  for  scatterer  characterization.  This  EM  field  analysis 
approach  that  is  essentially  based  on  scattering  vector  analysis  (rather  than  the  scattering  operator 
or  matrix)  is  similar  to  the  emerging  “compact  polarimetry”  methodology  [9] -[13]  that  exploits 
scattering  system  response  to  an  incomplete  set  of  input  EM  field  components.  Although,  the 
unique  scattering  transfer  function  is  ill-defined  mathematically,  certain  assumptions  may  be 
made  (e.g.,  symmetry)  to  construct  the  scattering  matrix  [  12]. 
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To  extract  the  most  information  about  the  target  from  its  random  backscatter  response  with 
unknown  orientation  relative  to  the  known  polarity  of  the  radar  illumination  (i.e.,  the  typical 
practical  scenario),  the  backscattering  measurement  potential  of  SAR  should  be  maximized.  The 
latter  is  maximized  if  and  only  if  the  data  products  are  the  4  Stokes  parameters  (or  their  logical 
equivalent)  [10].  This  is  in  agreement  with  the  statement  made  earlier  based  on  classical  EM 
physics  that  the  analytic  polarimetric  EM  signal  (polychromatic  and  partially  polarized)  is  best 
characterized  by  virtue  of  the  Stokes  vector  and  parameters.  The  above  statement  is  valid  for  a 
polarimetric  EM  signal  throughout  the  whole  spectrum,  from  the  microwave  to  visible  range. 

It  follows  from  the  above  discussion  that,  for  any  (or  regardless  of)  illumination  or  scattering 
source,  the  scattered  field  (and  hence  target)  characterization  is  optimized  by  using  the  4 -vector 
Stokes  parameters.  Hereby,  a  robust  platform  with  a  well-defined  set  of  observables  is  established 
for  interworking  between  EM  polarimetric  sensors  operating  at  different  regions  or  bands  of  the 
EM  spectrum,  e.g.,  SAR  and  optical  polarimetry. 

In  this  work,  the  wide -spectrum  scattered  EM  field  (partially/completely  polarized  or  coherent)  is 
formulated  in  a  unified  manner,  and  in  terms  of  components/products  of  significance  for 
scattering  analysis.  Since  these  EM  observables  are  characterizing  field  and  scattering  in  different 
regions  of  the  spectrum,  fusion  of  various  sensors  utilizing  EM  remote  sensing  for  advanced 
exploitation  is  envisioned. 
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3  Electromagnetic  (EM)  field  formulation  and 
observables 


3.1  EM  wave  analytic  description 

A  general  polychromatic  EM  wave  propagating  in  the  direction  r  can  be  analytically  represented 

by: 


£  (r,  0  =  E(r,  t)  exp  ( j(ut  -  k0L(  r,  t))) 


(2) 


E(r,f)  is  the  complex  vector  field  phasor, 


k0  — 


lit 


UJ 


where 

A  c 

wavenumber  and  speed  of  light  with  A  being  the  wavelength,  and 
position  vector.  Providing  that  the  waveform  bandwidth  is  small 
frequency,  a  quasi-monochromatic  model  can  be  adopted  as: 


and  c  are  the  vacuum 

r  —  (x,y,z)  denotes  the 
compared  to  the  center 


£  (r,  t)  =  E(f)  exp  ( j(u!t  -  0(f))) 


(3) 


In  the  spectral  domain,  the  polychromatic  EM  signal  (2)  can  be  represented  as: 


£  (r,  c o)  =  E(r,cu)  exp  (-j(k0L(r,  w))) 
= E  (r ,  'jj )  exp  (— j  (k  (r ,  w))) 


(4) 


In  (2)  and  (4),  L(r,t)  and  L(r,  u>)  represent  the  time  and  spectral  domain  eikonal  functions, 
respectively.  The  eikonal  function,  which  is  commonly  used  in  ray  theory,  characterizes  the  ray 
directions  and  the  wave  fronts.  The  time  invariant  (or  implicit  time  variant)  complex  vector  EM 
field/wave  in  the  EH  plane  (i.e.,  perpendicular  to  the  direction  of  wave  propagation  or  in 
direction  of  the  Poynting  vector)  may  be  expressed  as: 
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E(r)  =  Ex( r)  x  +  E  ( r)  y  -  Ex( r)  (x  +  p(r)  y) 


(5) 


where: 


P(r) 


Ey(r) 

E(r) 


(6) 


is  a  complex  factor  representing  the  polarization  ratio.  In  (5),  the  direction  of  wave  propagation  is 
assumed  to  be  in  the  z  direction  and  introduces  no  loss  of  generality.  Respectively, 

x,y  ( h , v  )  represent  the  scattering  or  polarimetric  basis  vectors.  One  should  note  that  the 

alignment  of  the  polarimetric  basis  vectors  x,  y  (  h,  v)  is  arbitrary.  The  x  alignment  may  be  set 
at  any  angle  a  (with  respect  to  the  universal  horizontal  plane)  and  the  y  alignment  can  be  found 

by  a  90-degree  right-hand  rotation  around  the  line  of  sight  (LOS),  i.e.,  a +  90.  It  is,  however, 
imperative  to  be  consistent  throughout  the  entire  formulation,  in  particular,  when  considering 
different  sensors. 

The  field  complex  vector  in  (5)  is  general.  It  can  be  partially  or  completely  polarized,  or 
completely  depolarized.  The  Stokes  vector  or  parameters  are  a  set  of  suitable  products  that  can 
comprehensively  describe  the  polarization  state  and  nature  of  EM  radiation  for  such  general 
waveforms.  The  Stokes  vector  associated  with  a  partially-polarized  EM  wave  (5)  is  given  by 

[1].[8],[14]: 


(e  *E  +E  E*) 

\  x  y  x  y  / 

i (e  *e  -e  e  *) 

J  \  X  y  X  y  / 


(7) 


where  (...)  represents  temporal  or  local  spatial  averaging.  The  first  vector  component  50  in  (7) 
represents  the  total  power  density  of  the  partially-polarized  field.  The  remaining  vector 
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components  5,  ,S7  ,S3  model  the  polarized  component  of  the  EM  field,  i.e.,  for  an  unpolarized 
EM  wave: 


sl=s7  =  s3  =  o 


(8) 


Decomposing  the  partially  polarized  EM  wave  into  completely  polarized  and  completely 
depolarized  field  components  using  (5),  gives: 


Ep(r)  =  £px(r)(x+Pp(r)y) 


(9) 


E(/(r)  =  JEc/x(r)(x+%(r)y) 


(10) 


so, 


ETotal  (r)=Ep  (r)+Eu  (r) 


(11) 


One  obtains  for  the  associated  Stokes  components: 


5 


u 

l 


=  0 


(12) 


and 


V  =2 


(Mr)|=1) 


(13) 
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Also, 


Alternatively, 


(15) 


One  may  also  parameterize  the  Stokes  vector  in  terms  of  the  EM  wave  polarization  ellipse  as: 
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*o' 

l 

Si 

—  J 

m  cos(2y)cos(2ri’) 

S2 

/ncos(2\)sin(20) 

SX 

m  si  n(2\) 

(16) 


where  la  denotes  the  total  power  of  the  partially -polarized  wave  (aligned  at  angle  a  as  described 
previously) ,  and  tn  is  the  degree  of  polarization  given  by  (using  (12)-(14)): 


s^-s"  ^i)2  +  fe)2  +  fe)2 

o  Total  o 


(17) 


In  ( 1 6),  x  and  'ip  are  the  characteristic  parameters  describing  the  EM  wave  polarization  ellipse, 
i.e.,  the  ellipticity  and  rotation  angles.  As  seen  from  (12)-(14)  and  (16)-(17),  the  characteristic 
parameters  %  and  if)  describe  the  polarized  portion  of  the  EM  wave.  As  such,  the  wave 
polarization  and  hence  the  Stokes  vector  can  be  mapped  on  to  a  3-dimensional  manifold  with 
space  angles  (x,V0  •  The  completely-polarized  waves  will  be  located  on  a  sphere  of  radius 
Ia  (i.e.,  the  Poincare  sphere),  while  the  and  partially  polarized  waves  lie  inside  the  sphere. 
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It  is  useful  to  have  the  analytic  representation  of  some  Stokes  vector  related  parameters  and 
products  for  a  partially  polarized  EM  waveform. 


The  degree  of  linear  polarization  is: 


m,  — 


J(S,?+(SJ 


—m  cos(2\) 


(18) 


The  circular  polarization  ratio  is: 


S0-S3J  m  L  tan(2;\) 

S0  +  S3  1  +  777  L  tan(2\) 


1+  a/h?2  —mL2 


(19) 
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The  polarization  ratio  phase  is: 


6  =  arc  tan 


(20) 


The  polarization  ratio  magnitude  angle  is: 


2  ap  =arccos 


Si  ' 
mS0/ 


(tan(ap)  =  |p|) 


(21) 


And: 


2x  =  arctan 


s3 

=  arctan 

mLS0 

(22) 


2-0  =  arctan 


(23) 


3.2  Stokes  vector  transformation 

As  addressed,  one  important  application  of  the  Stokes  EM  field  representation  is  to  study  and 
characterize  scattering  phenomena  as  a  result  of  target  interaction  with  the  EM  wave.  Therefore, 
analytic  transformation  of  the  Stokes  vector  associated  with  an  EM  waveform  transform  under 
certain  operations  (e.g.,  a  scattering  operator,  or  a  basis  transformation)  is  invaluable.  As  will  be 
more  evident  later,  the  Stokes  vector  and  parameter  transformations  play  an  important  role  in 
unifying  and  fusing  scattered  EM  wave  (partially  polarized)  characterization  in  different  spectrum 
bands  and  regimes. 


10 
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For  an  EM  polychromatic  wave  undergoing  a  transformation: 


E,(r)  =  f(r)©E,(r) 


(24) 


where  T  represents  the  operator  (e.g.,  Dyadic)  and  the  symbol  “0  ”  denotes  the  generalized  dot 
product  operator.  One  can  derive  the  Stokes  vector  transformation  operator  (using  straightforward 
operator  algebra)  as: 


G(r)=U,[fa(r)f(r))ur 


(25) 


where  Ul ,  U r  are  the  operators  defining  the  transformation  from  field  coherency  to  Stokes 
vector,  and  T“{r)  is  the  adj oint  operator  to  7"(r) .  In  a  more  familiar  matrix  format: 


^(r)(4x4)  U 


/(  4x4) 


r2x2  (r)  ®  T7 


2x2 


u 


(4x4) 


r(4x4) 


(26) 


In  (26),  the  symbol  0  denotes  the  Kronecker  matrix  product  and: 


U 


r 


10  0  1 
10  0  1 
0  110 
0  i  -i  0 


(27) 
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( for  general  transfromation ) 


U, 


10  0  1 
10  0  1 
0  110 
0  j  -j  o 


(28) 


1  0 
1  0 
0  1 
0  —j 


0  1 
0  1 
1  0 
j  0 


( for  backscattering  reception ) 


(29) 


Using  (25)-(29),  the  Stokes  matrix  transform  becomes: 


S,(r)=G(r)©Sf(r) 


(30) 


Upon  transformation  of  the  input  Stokes  vector,  various  alignments  of  the  scattered  field  Stokes 
vector  can  be  extracted  and  analyzed.  Projection  of  the  Stokes  vector  into  a  certain  direction  or 
alignment  is  not  straightforward  due  to  the  cross  coupling  of  the  field  components  (i.e.,  cross 
terms  in  power).  However,  through  definition  of  cross  Stokes  vectors  and  Stokes  basis  vectors 
(developed  and  to  be  reported),  projected  Stokes  and  power  terms  for  an  arbitrary  alignment  can 
be  determined. 


3.3  Examples 

Consider  the  well-known  example  of  a  trihedral  comer  reflector.  The  scattering  matrix  is  defined 
by: 


s  = 


1 

0 


0 

1 


(31) 
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According  to  (26),  the  Stokes  vector  transformation  matrix  is: 


10  0  0 

G=°  1  °  ° 

0  0  10 

0  0  0  -1 


(32) 


For  a  general  input  polarization,  the  Stokes  vector  is  parameterized  by  (see  (16)): 


S;  = 


1 

cos(2x,-)cos(2Vi) 

cos(2y.)sin(2^) 

sin(2x,) 


(33) 


and  the  scattered  Stokes  vector  is  readily  found  as: 


1 

1 

c°s(2xs,)cos(2'0s) 

cos(2x;-)cos(2^) 

cos(2xs,)sin(2?4) 

cos(2xI)sin(2'0I) 

sin(2  xs) 

-sin(2x,.) 

Projecting  this  output  vector  into  the  2  co-pol  and  cross-pol  directions  (using  the  procedure 
addressed  above),  one  obtains: 


Y  —  - 

co—  pol 


1+  cos(4%) 


(35) 
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and 


'  cross— pol 


1  —  cos(4\) 
2 


(36) 


The  incident  field  index  “i”  has  been  dropped  in  (35)  and  (36).  Figures  2  and  3  depict  the  co-pol 
and  cross -pol  response  to  a  general  polarization  characterized  by  ellipticity  and  rotation 
angles  (  tp)  .  This  behaviour  of  the  scattered  wave  from  the  trihedral  scatterer  can  be  compared 
to  its  scattering  characteristics  in  another  domain  such  as  EO  (by  using  the  described  EM  analytic 
framework)  for  advanced  exploitation. 


Figure  2.  Co-Pol  Reflection  ratio 
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Cross-Pol  Trihedral  Response 


Figure  3.  Cross-Pol  Reflection  ratio 
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4  Compact  Polarimetry  Concept 


The  identity  of  a  target  interacting  with  the  EM  field  and/or  wave  is  established  via  formulation 
of  scattering  and  by  virtue  of  a  scattering  operator  or  matrix.  Although  as  theoretically  discussed 
earlier,  performance  of  such  fully-polarized  radar  system  (i.e.,  quadrature-polarized)  is  unique 
and  has  no  substitute  in  providing  complete  target  backscattering  information,  it  comes  with  an 
attendant  cost,  e.g.,  lower  radar  swath  coverage,  higher  antenna  transmitter  power  requirements. 
As  a  result,  the  concept  of  “partial  polarimetry”  or  “compact  polarimetry”  has  emerged  in  recent 
years  and  several  papers  on  the  subject  have  been  published  [10]-[12],  [15]-[16],  The  main 
objective  of  this  line  of  research  is  to  achieve  certain  appealing  characteristics  of  a  fully-polarized 
system  without  actual  realization  of  a  quadrature -polarized  system.  Partial  or  dual  polarimetry  is 
essentially  a  step  up  from  a  single  channel  system  towards  the  fully  polarimetric  system  that  is  an 
effective  strategy  when  the  polarimetric  system  resources  are  limited  or  not  available.  It  is  also 
compatible  as  an  optional  mode  for  a  fully  polarimetric  radar  system  [9]. 

As  addressed,  the  backscattering  measurements  using  polarimetric  sensors  operating  at  any 
wavelength  or  frequency  are  optimized  and  the  extracted  information  maximized  using  Stokes  or 
equivalent  to  Stokes  parameters.  It  has  been  demonstrated  that  partial  or  compact  polarimetry  can 
emulate  many  aspects  of  the  full  or  quadrature -polarized  system  output  products  and  information 
[12], [16],  e.g.,  scattering  or  coherency  matrix,  by  incorporating  certain  a  priori  information  such 
as  reflection  and  rotation  symmetry.  Nonetheless,  Stokes  vector  parameterization  and  analysis 
represent  the  identifying  nature  and  characteristics  of  compact  or  dual-pol  polarimetry. 

In  a  compact  polarimetry  scenario  or  dual  polarization  radar/SAR  system,  a  single  polarization 
EM  signal  is  transmitted  (illumination)  and  the  backscattered  EM  signal  is  coherently  received  in 
two  orthogonal  polarizations.  The  choice  of  the  transmit  signal  is  arbitrary,  i.e.,  any  polarization, 
as  well  as  the  choice  for  two  orthogonal  receive  basis  polarizations  that  construct  a  complete  2- 
dimensional  space  to  characterize  the  scattered  (in  general,  partially  polarized)  EM  signal. 

Judicious  selection  of  these  polarization  sets  for  transmission  and  reception  can  optimize  the 
system  architecture  and  satisfy  the  operational  requirements  for  applications  of  interest.  Any 
assumed  orthogonal  polarization  set  at  the  receiving  end,  maximizes  the  backscattering 
information  through  constitution  of  a  complete  vector  space.  However,  choice  of  polarization 
basis  at  reception  can  affect  and  optimize  a  radar’s  design  with  respect  to  reliability,  architecture, 
mass  and  power  considerations.  The  nature  of  the  application  and  targets  of  interest  determines 
the  optimal  transmitted  polarization.  For  instance,  applications  manifesting  rotational  invariance 
imply  a  requirement  for  transmitted  circular  polarization  (CT).  Circular  polarization  transmit  and 
linear  (i.e.,  h  and  v)  polarization  receive  mode  (CTLR)  has  been  found  to  be  the  optimum 
architecture  for  lunar  or  planetary  exploration  and  a  good  alternative  to  linear  polarization  for 
Earth- observing  SAR  systems  [9].  One  of  the  initial  compact  polarimetric  modes  referred  to  as 

7T 

the  —  mode  [11],  [17],  utilizes  a  transmitted  field  polarized  at  45°  and  linear  polarization 

7T 

receive  orthogonal  set  (e.g.,  h,  v).  The  —  mode  is  successful  in  decomposition  analysis 
involving  scenes  with  scatterers  predominantly  oriented  in  horizontal  and  vertical  directions.  Both 
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of  these  mixed-polarity  modes,  CTLR  and  —  ,  enjoy  extensive  meteorological  radar  heritage 

4 

[  1 8] -[1 9] .  Also,  end-to-end  circular  polarimetry,  i.e.,  Circular  polarization  transmit  and  Circular 
polarization  receive  mode  (i.e.,  dual  circular  polarimetric  (DCP))  has  an  extensive  heritage  in 
radar  astronomy  [20]- [23], 


The  scattering  vectors  associated  with  the  described  modes  CTLR,  DCP,  and  — 

4 


expressed  in  terms  of  the  standard  quad-pol  scattering  matrix  components  as  following: 


can  be 


^■ctlr  — 


1 


^fih  JShv 

jSvv  +  s 

vh 


(37) 


k 


DCP  ~ 


$  RR 

Srl 


(38) 


where  “R”  and  “L”  denote  right  and  left  circular  polarizations,  and: 


k 


7 r 


4 


J_  shh  +  s 
svv  +  s 


hv 

vh 


(39) 


The  described  compact  polarimetry  fundamentals  are  similar  to  the  theoretical  foundation  that  the 
present  mutisensor  interworking  and  fusion  proposal  is  based  on.  In  brief,  the  core  methodology 
is  to  enhance  the  extracted  backscattering  information  by:  1 )  modeling  and  full  characterization 
of  the  scattered  EM  field  using  two  coherent  orthogonal  polarizations,  chosen  based  on  the 
system  design  and  architecture  requirements;  and  2)  designing  a  compact  polarimetry  mode  for 
intended  applications  based  on  respective  requirements  and  constraints.  In  the  current  report,  the 
described  methodology  and/or  analytic  tool  is  intended  to  formulate  and  enhance  the  common 
information  in  the  polarimetric  backscattered  field  from  different  sensors  (c.g..  radar/SAR  and 
electro-optical)  for  advanced  interworking  and  fusion. 
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5  EM  field  intensity  products  and  observables 


It  can  be  shown  that  the  Stokes  vector  associated  with  a  scattered  EM  wave  (16)  is  constructed  by 
measuring  the  field  intensity  at  4  polarizations  or  alignments  7a  ,1^,1  ,Ia  .  The  alignment 

angles  a,/3,ry,a  may  conventionally  be  set  at  0°,45°  ,90°, 135°  for  ease  of  computations. 
The  choice  of  these  4  angles  is  arbitrary  but  may  involve  more  complicated  computations. 

One  can  show  that  under  any  polarization  rotation,  the  first  Stokes  parameter  S0  is  invariant  and 
equal  to  sum  of  intensities  of  2  orthogonal  alignments: 


S0=Ie+Ie+9o  ( for  all  9) 


(40) 


Thus: 


S0  =  : 


L+I, 


a +90 


+  +  +  ++90 


(41) 


or  in  the  special  case: 


So- 2  (A) +  *45  +^90  +  J 135  ) 


(42) 


Other  Stokes  parameters  can  be  derived  using  a  similar  approach.  More  specifically: 


Si  —  Jq  J90 


(43) 


+  _  J45 ' 


135 


(44) 
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and: 


S3  = 


1 

2 


(m  4^/0  +/45  +/9Q  +/135  j  +  2^m  +  4)  (/0  /90  +/45  /135) 
+2  m  (/0  /45  +  70  /135  +  /45  /90  +  /90  /135 ) 


I2 


(45) 


where  m  is  the  degree  of  polarization.  In  order  to  provide  a  complete  picture  of  the  scattered  EM 
field,  i.e.,  Stokes  4-parameters,  the  degree  of  polarization  m  in  (45)  must  be  evaluated.  This 
parameter  may  be  estimated  using  a  number  of  techniques  such  as  entropy-based  extraction  [17], 
or  measurement  of  the  polarization  ratio  angle  (21). 

The  described  Stokes  vector,  parameters  and  related  quantities  constitute  a  unified  frame  of 
observables  for  the  scattered  wave  (fully  or  partially  polarized)  analysis  throughout  a  wide  EM 
spectrum.  As  such,  it  can  be  used  to  formulate  the  interworking  and  fusion  of  SAR  and  EO 
polarimetry  for  advanced  target  exploitation. 

As  an  example,  for  the  observed  scattered  field  intensity  that  can  be  compared  and  fused  for 
enhanced  exploitation,  consider  the  incident  EM  field  with  a  general  polarization  described  by 
(33).  As  derived,  the  scattered  field  by  a  trihedral  comer  backscatterer  received  at  the  receive 
antenna  is  given  by  (34). 

From  (40)-(44),  one  can  write: 


h 


Sq+S, 

2 


(46) 


ho~ 


(47) 


M5  ' 


^0+^2 


(48) 


DRDC  Ottawa  TM  2008-217 


19 


A  35 


So-S2 


(49) 


2 


Using  (34),  one  obtains: 


1+  cos(2x)  cos(2i/’) 


(50) 


T  _  1  —  cos(2x)  cos(2i/>) 

790  -  — 


(51) 


hs  — 


1+  cos(2x)sin(2'i/’) 
2 


(52) 


I135 


1  — cos(2x)sin(2'0) 
2 


(53) 


Simulations  depicted  in  Figures  4-7  show  the  above  scattered  field  intensity  variations  in  terms  of 
the  input  EM  field  polarization  variations  or  choice  of  transmitted  polarization,  i.e.,  the  ellipticity 
and  rotation  angles  ;y,U  . 

The  results  provide  a  tool  to  compare  or  combine  polarimetric  SAR  scattering  products  with  other 
EM  polarimetric  sensor  (e.g.,  EO)  products  based  on  a  robust  and  accurate  theoretical  foundation. 
The  significance  of  these  results  are  better  appreciated  when  considering  that  standard  EO 
polarimetric  products  are  typically  expressed  in  terms  of  4  field -intensity  outputs  (i.e., 
70  ,/45  ,/90  ,/135  )  of  polarization  filters  associated  with  the  described  4 -angles.  Hence,  the 
means  for  effective  fusion  is  provided. 
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Figure  4.  Trihedral  Corner  Reflector  Scattered  Field  Intensity  I0  (General  Illumination 
Polarization  \ f  ). 
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Figure  5.  Trihedral  Comer  Reflector  Scattered  Field  Intensity  /45  (General  Illumination 
Polarization  ). 
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Figure  6.  Trihedral  Comer  Reflector  Scattered  Field  Intensity  /90  (General  Illumination 
Polarization  x ,  f  ). 
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Normalized 

Intensity 


Figure  7.  Trihedral  Corner  Reflector  Scattered  Field  Intensity  /135  (General  Illumination 
Polarization  XFP  )■ 


24 


DRDC  Ottawa  TM  2008-217 


Normaliyd 
Intensity 


Figure  8.  Trihedral  Corner  Reflector  Scattered  Field  Total  Intensity  in  4  Manifolds,  i.  e. , 
I(j ,  /45  ,  /90  ,  /|  35  (General  Illumination  Polarization  XdlJ  )■ 
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6  Conclusions  and  future  direction 


Polarimetric  characteristics  of  the  EM  waveform  provide  valuable  additional  information  for 
target  recognition  compared  to  conventional  remote  sensing  techniques.  The  objective  here  has 
been  to  examine  the  fusion  of  EM  polarimetric  quantities  and  measurables  associated  with  certain 
classes  of  targets  in  different  regions  of  the  EM  spectrum.  One  prominent  example  is  the 
combination  of  radar/SAR  polarimetry  with  EM  polarimetric  sensors  in  the  electro- optical 
domain. 

The  described  EM  polarimetry  fusion  within  different  bands  must  be  carried  out  within  a  robust 
theoretical  framework  based  on  the  physics  of  partially  polarized  EM  waves  scattered  at  different 
operating  wavelengths.  Using  this  framework,  sensor  fusion  can  yield  a  significant  advancement 
in  target  recognition  due  to  the  following:  target  scattering  characteristics  (i.e.,  identity)  are 
better -established  in  a  wider  range  of  the  EM  spectrum,  even  if  this  spectrum  expansion  is 
incoherent  (i.e.,  incoherent  bandwidths);  polarimetric  backscatter  characteristics  for  certain  target 
classes  of  interest  are  consistent  throughout  various  regions  of  the  EM  spectrum  (e.g., 
polarization  selectivity  for  manmade  objects).  Fusion  of  EM  polarimetric  returns  in  different 
bands  can,  in  principle,  enhance  the  contrast  of  such  objects/targets  compared  to  those  with  low 
or  without  the  addressed  polarization  selectivity  (e.g.,  clutter).  The  described  EM  polarimetric 
nature  (i.e.,  polarization  selectivity)  manifested  in  different  regions  of  the  spectrum  can  be 
effectively  exploited  for  modification  and  calibration  puiposes. 

As  addressed,  the  analytic  framework  should  be  based  on  a  careful  study  and  modeling  of  the  EM 
wave  scattering  and  respective  characteristic  observables  in  different  regions  of  the  spectrum. 
C  onsidering  the  general  property  of  the  scattered  EM  wave  being  polychromatic  and  partially 
polarized,  the  analytic  polarimetric  EM  signal  can  be  characterized  by  the  virtue  of  the  Stokes 
vector  and  parameters.  These  parameters  and  related  products  can  be  used  as  an  effective  means 
for  sensor  data  combination  and  fusion. 

The  proven  potent  and  fast  growing  field  of  partial  or  compact  polarimetry  is  established  b  ased  on 
fundamentals  similar  to  the  theoretical  foundation  that  the  present  mutisensor  interworking  and 
fusion  proposal  is  based  on.  The  core  methodology  is  to  enhance  the  extracted  backscattering 
information  by: 

I.  modeling  and  full  characterization  of  the  scattered  EM  field  using  two  coherent  orthogonal 
polarizations,  chosen  based  on  the  system  design  and  architecture  requirements; 

II.  designing  compact  polarimetry  modes,  or  essentially  transmit  or  illumination  optimization 
for  intended  applications  based  on  respective  requirements  and  constraints. 

One  should  note,  however,  to  maximize  the  extracted  common  information  for  sensor  fusion,  full 
intensity  or  Stokes  parameters  (4-vector)  description  is  required  (as  addressed  in  Section  5). 

Here,  the  described  methodology  and/or  analytic  tool  is  used  to  formulate  and  maximize  common 
information  in  the  polarimetric  backscattered  field  from  different  sensors  (e.g.,  radar/SAR  and 
optical)  for  advanced  fusion. 
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To  provide  a  means  for  advanced  exploitation  and  recognition  through  an  optimal  combination  of 
sensor  data  (e.g.,  SAR  and  EO),  the  following  directions  are  proposed: 

1.  Secure  polarimetric  data  for  a  certain  target  class  (or  classes)  of  interest  using  radar  and 
EO  sensors; 

2.  Formulate/format  the  data  in  terms  of  analytic  EM  wave  vectors  (Stokes  vectors  and 
parameters);  and 

3.  Investigate  scattered  EM  wave  polarization  characteristics  related  to  the  physics  of  the 
target  that  are  consistent  in  all  (or  both)  EM  regions  for  enhanced  fusion  calibration  and 
exploitation  in  one  domain  using  the  data  and  target  characteristics  from  the  other 
domain. 
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